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I. Introduction 


As part of a project to produce a multi-year assimilated dataset, the 
Data Assimilation Office (DAO) in the Goddard Laboratory for 
Atmospheres has made extensive modifications to the existing general 
circulation model (the Goddard Earth Observing System Model, GEOS-1). 
An important aspect of the modifications was the development of improved 
lower boundary conditions. In particular, GEOS-1 does not have a coupled 
land surface model and requires soil moisture to be specified. The soil 
moisture, however, is not a quantity readily derivable from conventional 
observations without adequate surface radiation measurements. 

Thomthwaite (1948) proposed a procedure to compute soil moisture 
from observed monthly mean precipitation and estimated potential 
evapotranspiration using surface air temperature. The procedure has been 
tested in many field studies (Mintz and Serafini, 1984 and references 
therein) and shown to give a reasonable approximation to local 
measurements of the soil moisture and evapotranspiration. Mintz and 
Serafini adopted this procedure to generate gridded climatological soil 
moisture distributions which have been used as a part of the climatological 
GCM boundary conditions. In this study, we have used the same 
procedure with some modification to generate soil moisture fields for each 
month of the period 1979 - 1989. 

Section II describes the surface station climatological dataset. 
Section III describes the procedure to obtain soil moisture estimates, and 
Section IV lists the global quantities included in this document. 


II. Data 

A copy of the world monthly surface station climatology dataset was 
obtained from the National Center for Atmospheric Research (NCAR). 
The dataset includes monthly mean sea level pressure, station pressure, 
surface air temperature and monthly total precipitation at surface stations 
from mid-1700’s up to 1989. The number of station reports varies from 
year to year, but is on the order of 2500 in more recent periods. Most of 
the data included in this dataset were obtained from the National Climatic 
Data Center (NCDC), Asheville, North Carolina. However, several 
individuals contributed towards upgrading and completing the dataset. 
These include Sharon Nicholson (Florida State University) for the African 
precipitation data and Dennis Shea (NCAR) for data obtained directly from 
various countries. A detailed description of the dataset and data sources 
can be found in Spangler and Jenne (1990). 
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For the present work, a subset of the monthly mean surface air 
temperature and the monthly total precipitation was created for the years 
1979 - 1989. For any year, stations with more than 6 months of missing 
data were not used. The number of stations included in the soil moisture 
computation is listed in Table 1. An example of data coverage is given in 
Fig. 1; the crosses indicate the location of stations with more than 6 
months of data in 1985. 


m. Soil Moisture Estimate 

Mintz and Serafini (1984) derived a procedure to compute soil 
moisture based on the method proposed by Thomthwaite (1948). An 
equation for the soil moisture, w, is written as, 

3w _ * ... 

— = P-E , Wmax =W , (1) 

dt 

where 

E = EpxP(w,w*) . (2) 


Here, P denotes the precipitation rate, E the rate of evapotranspiration, w 
the maximum soil moisture available for evapotranspiration which is defined 
as the difference between the moisture storage capacity of soil and soil 

moisture at the wilting point of plants, Ep the potential evapotranspiration, 

and P the evapotranspiration coefficient. The magnitude of w is fixed at 

150 mm as in Mintz and Serafini, and any amount greater than w is 
treated as runoff. 

Following Thomthwaite (1948), the potential evapotranspiration for a 
month is estimated as a function of the monthly mean surface air 
temperature and duration of daylight. 


E P = 


0 , 

16L(10T a / I) a , 

-415.85 + 32.35T a - 0.43T A 2 , 


T a < 0°c 

0 < T A < 26.5° c (3) 
T a > 26.5° c 


where 


I = I (Tab, / 5) 1514 , 

m=l 
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a = (6.75xl0‘ 7 I 3 ) - (7.71xl0* 5 I 2 ) + (1.79xl0’ 2 I) + 0.49 , 
and 

L = (D/30)(h/12) . 

Here D is the number of days in the month, h the number of hours of 

daylight, and Ta the monthly mean surface air temperature in each month 
of the year. 

The evapotranspiration coefficient, P , is given as 

P = 1 - exp[-6.8(w/ w*)] . (4) 


The above relationship is based on the analysis of field measurements by 
Nappo (1975) as adopted by Mintz and Serafim (1984), and represents the 

curve of best fit to (E/Ep) plotted against the volumetric soil moisture 
content. 

Although the above procedure to estimate the evapotranspiration is 
somewhat cmde, it is shown to give reasonable results and has been used 
extensively (Mintz and Serafini, 1984 and references therein). In some 
cases, however, it is known to strongly underestimate the potential 
evapotranspiration (Sellers, 1965). 

Mintz and Serafini (1984) developed an iterative procedure to solve 
(1) for w and E based on climatological monthly mean precipitation and 
surface air temperature. In the present computation, an analytic solution 
was obtained after substituting (2) and (4) into (1). The solution is given as: 


w _ _ (P-Ep)exp{-q(P-E p )t) 

a (P - E p )exp(aw 0 ) + E p [l-exp{-a(P-E p )t}] 

where 

a = 6.8 / w* 


(5) 


Here, wo is the initial soil moisture and time t is elapsed time. For each 

month, the soil moisture, w, was computed daily with an averaged 

precipitation rate and potential evapotranspiration based on the monthly 

means. The monthly mean soil moisture was obtained by averaging the 

daily values for the month and the evapotranspiration E was estimated from 

(2). As in Mintz and Serafini, runoff was accounted for by not allowing w 
* 

to exceed w . 

It should be noted here that the above procedure does not account 
for frozen ground and ice formation. Since there is no evapotranspiration 



when the air temperature is below the freezing point, the precipitation is 
simply accumulated as soil moisture. This might overestimate the soil 
moisture during the winter season. 

The soil moisture was computed at each station for the period of 
1979 - 1989 in one year segments. For each year, the solution was allowed 
two months of ’spin-up' period by starting the integration at the beginning of 
November of the previous year as initial conditions. For convenience, 
November climatological soil moisture values of Mintz and Serafim were 
used as the initial conditions, since several experiments with the initial 
conditions indicated that the results were independent of the initial 
conditions after two months of integration. The procedure was also 
reinitialized with climatological soil moisture values after periods of missing 
data. 

In order to generate data suitable for GCM boundary conditions, the 
station values of soil moisture were interpolated to a 2x2.5 degree 
latitude/longitude grid. The interpolation was performed by averaging 
station values within a 300 km radius of each grid point. The averaging 
weights are proportional to the inverse of the square of the distance of 
stations from the grid point. The value at a grid point with no station data 
within the 300 km radius was denoted as undefined. 


IV. List of Quantities 

The following is a list of global maps included in this report. 

1) Monthly mean ground wetness defined as the ratio 
of soil moisture to the maximum available moisture 

(w/w*, w* = 150 mm) for the period 1979 - 1989. 

2) 11 -year average of monthly mean ground wetness, 
precipitation, surface air temperature, evapotranspiration 
and potential evapotranspiraton. Also included is the 
standard deviation of the monthly mean ground wetness 
based on the 11 years of data. 

For the plots, precipitation, surface air temperature, 
evapotranspiration and potential evapotranspiration were interpolated in 
the same manner as described in Section in. On the maps, areas with 
undefined values are indicated by hatching and the interpretation for 
different shadings is included where it is needed. For the standard 
deviation of ground wetness, areas with values greater than 0.05 are lightly 
shaded. 

The datasets described in this document and the gridded monthly 
mean surface air temperature, monthly total precipitation, potential 
evapotranspiration and evapotranspiration for the period 1979 - 1989 are 
available to the general community. For further information on these 
datasets, contact 

Head, Data Assimilation Office 

Code 910.3 

NASA/Goddard Space Flight Center 

Greenbelt, MD 20771 . 
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Table 1. Number of surface stations included for soil moisture estimates 


Year No. of stations 


1979 

1351 

1980 

1345 

1981 

1261 

1982 

1215 

1983 

1320 

1984 

1254 

1985 

1151 

1986 

1127 

1987 

1369 

1988 

1426 

1989 

1343 


IX 



STATION LOCATION FOR SURFACE DATA IN 1985 



X 


LONGITUDE 

Fig. 1 : An example of the data coverage for soil moisture calculation. The crossi 
indicate the location of stations with surface data for more than 6 months in 1985. 
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